Stably stratified magnetized stars in general relativity 
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We construct magnetized stars composed of a fluid stably stratified by entropy gradients in the 
framework of general relativity, assuming ideal magnetohydrodynamics and employing a barotropic 
equation of state. We first revisit basic equations for describing stably-stratified stationary ax- 
isymmetric stars containing both poloidal and toroidal magnetic fields. As sample models, the 
magnetized stars considered by Ioka and Sasaki [231] - inside which the magnetic fields are confined, 
are modified to the ones stably stratified. The magnetized stars newly constructed in this study 
are believed to be more stable than the existing relativistic models because they have both poloidal 
and toroidal magnetic fields with comparable strength, and magnetic buoyancy instabilities near 
the surface of the star, which can be stabilized by the stratification, are suppressed. 

PACS numbers: 04.40. Dg,97.60.Jd 



I. INTRODUCTION 

Recent observations established that soft-gamma re- 
peaters (SGRs) and anomalous x-ray pulsars (AXPs) 
are the so-called magnetars, i.e., highly magnetized neu- 
tron stars whose surface field strength is as large as 
~ 10 14 — 10 15 G [ij-il- The presence of the magnetars has 
activated studies on equilibrium configurations of mag- 
netized stars, which have a long history. 

Since the pioneering work by Chandrasekhar and 
Fermi [f| , an enormous number of studies has been done 
for exploring structures of magnetized stars: Prcndcr- 
gast 0] and Woltjer [H calculated equilibrium configura- 
tions of magnetized stars having mixed poloidal-toroidal 
fields, where the magnetic fields are treated as first- 
order perturbations around a spherical star (see, also, 
Ref. [9(). Monaghan [lOj studied magnetized stars con- 
taining purely poloidal magnetic fields. Ioka [ll| devel- 
oped the works by Prendergast and Woltjer into those 
of second-order perturbations to study magnetic effects 
on the stellar structures. He also employed the re- 
sults obtained to explain magnetar activities. Miketinac 
obtained magnetized stars containing purely toroidal 
fields [l2j and purely poloidal fields [l3[ by solving exact 
master equations numerically. Tomimura and Eriguchi 
developed a numerical method for obtaining magnetized 
stars with mixed poloidal-toroidal fields usin g a non- 
perturbative technique [3 (see, also, Refs. [IB-[I3l). 
Duez and Mathis variationally considered the lowest- 
energy equilibrium states for a fixed magnetic helicity 
and constructed equilibria of magnetized stars having 
mixed poloidal-toroidal fields by a perturbation tech- 
nique [18j |. 

General relativistic models of magnetized neutron stars 
have been also explored. Bocquet et al. and Cardall 
et al. 20] obtained relativistic neutron star models with 
purely poloidal magnetic fields. Using a perturbativc 
technique, Konno et al. [2lJ calculated similar models. 
Kiuchi and Yoshida [22| computed magnetized stars with 
purely toroidal fields. Ioka and Sasaki (23|, Colaiuda et 
al. 124], and Ciolfi et al. 125, 26] derived relativistic stellar 



models having both toroidal and poloidal magnetic fields 
with a perturbative technique. Although progress has 
been achieved in this field, further studies are required 
because all the magnetized star models are constructed 
by some special magnetic-field configurations which may 
not be realistic. In particular, it is not clear at all whether 
their models are stable. 

The stability of the magnetized star is an important 
issue, because only stable equilibrium models are viable. 
Stability analyses of magnetized stars have been per- 
formed by many works, since the pioneering work by 
Tayler [27| , who showed that stars having purely toroidal 
magnetic fields are unstable. Wright [28| subsequently 
showed that there is the same type of the instability, 
the so-called pinch-type instability, for stars containing 
purely poloidal magnetic fields. He also suggested the 
possibility that stars having mixed poloidal-toroidal mag- 
netic fields may be stable if the strength of both com- 
ponents is comparable (see, also, Refs. [2!| Ass- 
che et al. [31| proved that the pinch-type instability in 
general arises in magnetized stars with purely poloidal 
fields, and Wright |28| and Markey and Tayler [29j stud- 
ied this instability for particular magnetic-field configu- 
rations. Flowers and Ruderman [32] found that another 
type of instability occurs in purely poloidal magnetic- 
field configurations. 

All those classical stability analyses have been done 
by a method of an energy principle in the framework of 
Newtonian dynamics (see, also, Refs. [H, H3])- Another 
approach is a local analysis, with which Acheson (35| in- 
vestigated the stability of rotating magnetized stars con- 
taining purely toroidal fields in detail in the framework 
of Newtonian dynamics (see, also, Refs. [36|, and 
derived detailed stability conditions for purely toroidal 
magnetic fields buried inside rotating stars with dissipa- 
tion. Note that although it is an approximate approach, 
the local analysis can take account of realistic effects on 
the stability like rotation, heat conduction, and resis- 
tivity, which cannot be included in a method of the en- 
ergy principle. Bonanno and Urpin analyzed the axisym- 
metric stability [38| and the non-axisymmetric stabil- 
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ity [391 ] of cylindrical equilibrium configurations possess- 
ing mixed poloidal-toroidal fields, while ignoring com- 
pressibility and stratification of the fluid. 

Recently the stability problem of the magnetized star 
has been approached from another direction. By fol- 
lowing the time evolution of small random initial mag- 
netic fields around a spherical star in the framework 
of Newtonian resistive magnetohydrodynamics, Braith- 
waite and Spruit [13, EH obtained stable configurations of 
a magnetized star that are formed as a self-organization 
phenomenon. The resulting stable magnetic fields have 
both poloidal and toroidal components with comparable 
strength and support the conjecture for stability condi- 
tions of the magnetized star given by the classical studies 
mentioned before. By using the numerical magnetohy- 
drodynamics simulation, Braithwaite 42] studied stabil- 
ity conditions for the magnetized stars and obtained a 
stability condition for his models given in terms of the 
ratio of the poloidal magnetic energy to the total mag- 
netic energy which is of order unity. Duez et al. showed 
that magnetized stars constructed in Ref. [l8| exhibit no 
instability for several Alfven time scales in their numeri- 
cal simulations [43| . This fact reconfirms the results given 
by Braithwaite. Lander and Jones explored the stability 
of magnetized stars by numerically solving the time evo- 
lution of linear perturbations around the stars in their 
series of papers [44M46j . For the purely toroidal/poloidal 
field cases, their results are consistent with those of the 
classical stability analysis, i.e., the pinch- type instability 
is observed near the symmetry and the magnetic axes for 
the purely toroidal and purely poloidal field cases, respec- 
tively. They also assessed the stability of various magne- 
tized stars with mixed poloidal-toroidal fields and found 
that all their models considered suffer from the pinch- 
type instability even for the cases in which the poloidal 
and toroidal components have comparable strength [46| . 
It is obvious that the results by Lander and Jones are 
incompatible with those by Braithwaite and his collabo- 
rators [iH HD]. Lander and Jones discussed the possibility 
that some physics missing in their study would suppress 
the instability they found. We infer that in particular, 
stratification of the fluid will be a key ingredient, which 
is taken into account in the analyses of Refs. [HI, [43} but 
not in the analyses of Ref. [46|. We will return to this 
point later. 

General relativistic magnetized stars have been also 
analyzed recently. By numerical-relativity simulations, 
Kiuchi et al. [43, [48| investigated the stability of the 
magnetized stars with purely toroidal magnetic fields 
obtained by Kiuchi and Yoshida [22[. They showed 
that the stars with some specific distributions of mag- 
netic fields are stable against axisymmetric perturbations 
but all the models considered are unstable against non- 
axisymmetric perturbations due to the strong magnetic 
buoyancy instability near the surface of the stars. The 
initial behavior of the instability observed in Refs. [471I48} 
is consistent with that expected by the Newtonian linear 
analyses by Acheson (35[. Lasky et al. and Ciolfi et 



al. [50l ] showed by numerical-relativity simulations that 
the purely poloidal magnetic fields, obtained by Bocquet 
et al. are unstable due to the pinch-type instability 
near the magnetic axis as predicted by the Newtonian 
linear analyses @, [29[ . 

All these recent general relativistic magnetohydrody- 
namics simulations have contributed a lot to the progress 
of the stability analyses of general relativistic magnetized 
stars. However, the numerical simulations have been per- 
formed for equilibrium stars composed of non- stratified 
fluids. The assumption of non-stratification is often used 
and quite reasonable as a first approximation for explor- 
ing cold neutron stars, even though the cold neutron star 
is expected to be highly stably stratified by the compo- 
sition gradient (see, e.g., Refs. (5lT[53l ]V If the effects of 
magnetic fields are taken into account for the neutron 
star models in their stability analysis, however, the situ- 
ation changes drastically. As discussed by many authors, 
e.g., Reisenegger [Hi} and Kiuchi et al. 48], it has long 
been known that the magnetic buoyancy makes the mag- 
netized star unstable and that the stable stratification is 
necessary to remove the magnetic buoyancy instability 
(the so-called Parker instability |55|). It should be em- 
phasized that in the assumption that the stellar matter is 
stably stratified, Braithwaite and Spruit [40|,|4l| obtained 
stable magnetized stars of simple magnetic configurations 
in their numerical simulations. Therefore, we infer that 
a stable stratification is one of the key ingredients for 
stable configurations of the magnetized stars. Note that 
Lander and Jones [46[ indicated the possibility that for 
stars containing mixed poloidal-toroidal magnetic fields, 
some weak instability associated with poloidal magnetic 
fields may not be removed by stable stratification (see, 
also, Ref. [38j). However, at the moment, no definite con- 
clusion has been obtained. A reason that non-stratified 
magnetized stars are employed for the stability analyses 
in general relativity is that no model of stratified magne- 
tized stars have been constructed, although in the frame- 
work of Newtonian dynamics, magnetized stars with sta- 
ble stratification due to composition gradients have been 
obtained [HfiUEl. 



In this paper, thus, we study stably stratified and mag- 
netized stars in the framework of general relativity aim- 
ing at giving a prescription for constructing them. First 
of all, we describe a general formulation to obtain sta- 
tionary axisymmetric magnetized stars composed of both 
toroidal and poloidal magnetic fields with stratification 
due to entropy gradients assuming ideal magnetohydro- 
dynamics and employing a barotropic equation of state. 
As sample models, the magnetized stars considered by 
Ioka and Sasaki [23| , which contain poloidal and toroidal 
fields of comparable strength only inside the stars, are 
modified to be the ones stably stratified. Note that to 
date, no general relativistic magnetized stars containing 
mixed poloidal-toroidal fields have been constructed with 
a non-perturbative approach because of difficulties in the 
treatment of non-circular spacetimes (see, e.g., Ref. 159} ). 
Finally, we describe the reason that the magnetized stars 
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obtained in the present study are more stable than the 
existing relativistic models. 



II. BASIC EQUATIONS FOR THE GENERAL 
RELATIVISTIC IDEAL 
MAGNETOHYDRODYNAMICS 

We consider perfect fluids coupled with electromag- 
netic fields, described by the basic equations summarized 
as follows. The baryon mass conservation equation: 



V„(ptt") = 0, 



(1) 



where p and are the rest-mass density and the fluid 
four velocity, respectively with V M being the covariant 
derivative associated with the metric g^. The two sets 
of the Maxwell equations: 



V„iW = 4ttJ' j , 



0. 



(2) 
(3) 



where F^ v and J M are the Faraday tensor and the current 
four vector, respectively. The total energy-momentum 
conservation law: 



V„T^' y = , 



(4) 



where T^ v is the total energy-momentum tensor, defined 
by 



T" v = phu^u" '+Pg^ ' + — 

47T 



F» a F" a - ^ v F al3 F aP 



(5) 

with h and P being the specific enthalpy and the pres- 
sure, respectively. Here, the specific enthalpy is written, 
in terms of the specific internal energy e, the pressure P, 
and the rest-mass density p, by 



h = 1 + e + P/p . 



(0) 



It is convenient to introduce the electric field and 
the magnetic field observed by an observer associated 
with the matter four velocity , defined by 



= F^u h 
1 



B„ 



(7) 
(8) 



where e M „ aj 3 is the Levi-Civita tensor with 60123 



-9 

with g being the determinant of the metric g^. In a 
large number of intriguing astrophysical problems, fluids 
coupled with electromagnetic fields can be approximated 
with perfectly conductive ones. Thus, we further assume 
the condition of perfect conductivity, 



E„ 



F^u" 



0. 



The dual tensor of F^" is then 



f 



p*^ = e^Fap = B»u v - B 



(9) 



(10) 



Equation §4$ is often decomposed into two sets of equa- 
tions, the energy equation and the momentum equation 
in the fluid rest frame, respectively, given by 

-u„V v T" v = u v V v {p{l + e)} + phV u u v 

= pu v V v e + PV v u v = , (11) 
g^VJ"" = phu v V uUtl + g£V„P - F^r = Qjl2) 

where q^ v = g^ v + u^u v . The perfect conductivity con- 
dition has been used to derive Equations (| 1 1 1) and (|12l) . 
Using Equation ([TJ and the first law of thermodynamics, 



P 

de = ^rdp + TdS. 
P 2 



(13) 



Equation pip is recasted into the entropy equation, 

u"V u S = 0, (14) 

where S and T are the specific entropy and the temper- 
ature, respectively. 

To construct a magnetized star, we employ barotropic 
equations of state given by 



P = P(p), e = e(p) 



(15) 



This equation of state is often used when studying a re- 
alistic star. 

It should be noted that non-stratified isentropic fluids, 
characterized by de = p~ 2 Pdp, have been assumed in a 
large number of studies of neutron stars in general rela- 
tivity. When studying cold neutron stars in a chemical 
equilibrium, this simplification may be accepted. How- 
ever, this is not the case if effects of magnetic fields are 
taken into account. As pointed out by many authors 
(see, e.g., Refs. |54( and [48j). the stable stratification is 
a necessary condition for magnetized stars to be stable. 
The reason is that magnetic flux tubes inside the star 
basically suffer from magnetic buoyancy and are forced 
to move toward their surface, resulting in a desperate in- 
stability of the stars. To counteract this magnetic buoy- 
ancy and to stabilize the stars, the stable stratification is 
required. For sufficiently cold neutron stars, the proton- 
neutron composition gradient is a candidate for the strat- 
ification [5l|. When studying on the structure of a mag- 
netized star, therefore, it is crucial to take into account 
this effect. Thus, the condition, de — p~ 2 Pdp, which im- 
plies that the star is non-stratified, should not be a priori 
assumed. Note that whether the star is barotropic or not 
is independent of whether the star is stably stratified or 
not. With the one-parameter equation of state (fT5|) . it is 
possible to have a stably stratified star. 



III. MASTER EQUATIONS FOR EQUILIBRIUM 
MAGNETIZED STARS WITH STABLE 
STRATIFICATIONS 

In this section, we derive the master equations for 
describing stably stratified axisymmetric rotating stars 
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composed of mixed poloidal-toroidal magnetic fields. We 
take the time and azimuthal coordinates as x° — t and 
x 3 = ip, respectively. Then, components of the two 
Killing vectors may be written as — Sq and ip^ = S 3 , 
where 5$ is the Kronecker delta. The other two spa- 
tial coordinate variables are written as x 1 and x 2 in this 
section. Thus, the quantities describing the equilibrium 
stars are basically functions of x 1 and x 2 only. Hence- 
forth, capital Latin indices (A, J5,C, •••) run from 1 to 
2. 

Because of the assumption of the axial symmetry and 
stationarity, Equation (|14[) is written as 



i 7 , 



12 



u A d A S = 0. 



(16) 



This equation means that the specific entropy has to be 
constant along streamlines on the stellar meridional plane 



unless 



0. However, the constant specific entropy 



distribution along streamlines is not realized for stable 
magnetized stars because of the following reason: For 
stationary axisymmctric stars, streamlines on the merid- 
ional plane, in general, are closed curves. Thus, it is 
inevitable that there exits an unstably stratified region 
(convectively unstable region) where (V /i P)(V M S < ) < 
is satisfied. To construct a stably stratified (convectively 
stable) star, thus, we have to assume 



u A = 0. 



(17) 



Then, we have no condition for S apart from the as- 
sumption of the stationarity and the axial symmetry. In 
other words, we can freely choose a functional form of S 
if u A — is assumed. 

The assumption, u A = 0, is an essential difference be- 
tween our study and the study of Ref. [23| in which the 
isentropic meridional flow is taken into account. The 
fluid four-velocity is, thus, given by 



7(f + Qip") 



(18) 



where O is the angular velocity of the fluid and 7 = u° . 

From Equations (O, and the integrability condi- 
tion for the momentum equation (|12[) . we have the fol- 
lowing relations: 



F A3 



d 

dx A 



Ft 



OA 



F 03 = : 



V^gF 12 = f (*) , 

Q F (^)u°+u 3 = 0, 



In 7 



dP 
ph 



/u(*)d* + C = . 



(19) 

(20) 
(21) 
(22) 
(23) 

(24) 



=?', (25) 



j 3 - n F j° = ph{fi(*) + ~fu 3 n' F } + 



where fl F , T, and pt are arbitrary functions of the flux 
function it, which is the azimuthal component of the vec- 
tor potential associated with F^ as shown in Equa- 
tion (I19|) . Here, C is an integral constant, and the prime 
denotes the derivative with respect to ^ . Note that Cl F 
is sometimes called the angular velocity of the magnetic 
field line. Substituting Equation (|18|) into Equation 
we obtain 



fl F = CI. 



(26) 



Equation (|24)) corresponds to the equation of the hydro- 
static equilibrium. 

For cases of non-stratified stars, the present formula- 
tion may be derived from that of Ioka and Sasaki 23] by 
taking a limit. First, we focus on Equations (16) - (18) 
of Ref. [23j . Using a relation 



7(u + Clu 3 ) 



1 



(27) 



or u^Ufj_ = —1, we may rewrite Equation (f2~4"|) as 

(u + flu 3 )eJ P h = -eJ , (28) 

which corresponds to Equation (16) of Ref. 23]. We 
therefore found that functions p and D in Ref. [23j j corre- 
spond toexp(J j^-) and exp(J /i(^')(i^ , — C), respectively. 
Note that p, of Ref. [23[ is exactly the same as h of the 
present study and that h~ 1 dh = (ph)~ 1 dP in the isen- 
tropic case de = p~ 2 Pdp. Equations (17) and (18) of 
Ref. [23| are, respectively, derived from the energy and 
angular momentum conservation equations, 



V M (TV) = 0, V M (T'W)=0. 
These may be explicitly written by 



0. 



(29) 



(30) 



with being or ^ . Since we assume the condition 
(fl"7f , the first term in the left-hand side of Equation (|30|) 
automatically vanishes and we obtain two relations, 



X 



Y 



(31) 



where X and Y are arbitrary functions of the flux func- 
tion "J. These relations were already derived in Equations 
([21]) and (J22]). Thus, we find that X = f fl and Y = f . 
The functions X and Y are, in terms of the functions of 
Ref. [1^1, C, E, and L, given by 



x =- 4 4> 

From Equation (20) of Ref. 
flL/C. As argued in Ref. 



Y = 



. L 
-An-. 



(32) 



23 

23] 



wehave£)/C = S/C- 
the no-meridional flow 
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limit of magnetized stars with mixed poloidal-toroidal 
magnetic fields is given by the limit C — > oo with D, L/C 
and E/C kept to be finite. Thus, we find that D/C = 
E/C-QL/C becomes = -4ttE/C+4ttQL/C = X-ClY 
in this limit, which is automatically satisfied by Equation 
(f3"Tj) in the present situation. 

Once the metric coefficients are given, all the compo- 
nents of are written as functions of ^ and / dx A 
through the two arbitrary functions of the flux func- 
tion, f2 F and T. Using the Maxwell equation ©, thus, 
J° — ftp J 3 appearing in the left-hand side of Equa- 
tion (|25l) can be written as a second-order elliptic-type 
partial differential operator for ^(rr 1 , x 2 ). Equation ([3]) 
in conjunction with Equation (|25|) may then be solved 
to obtain a distribution of the magnetic fields around a 
star. This equation is often called the Grad-Shafranov 
(GS) equation [&(| when it is written as a partial differ- 
ential equation for vp. 

It is useful to introduce some global quantities to char- 
acterize equilibrium solutions of stars. For equilibrium 
states of magnetized stars, the total baryon rest mass 
M* , the internal thermal energy E lnt , and the electro- 
magnetic energy .Eem may be defined as 



M* = 



-gd x, 



Eem = ^ / B»B^^g~d s x 
on 



(33) 
(34) 
(35) 



(see, e.g., Ref. [22 



IV. MAGNETIC FIELDS AROUND A 
SPHERICAL STAR AND THEIR EFFECTS ON 
THE STELLAR STRUCTURES 

A. Spherical stars with no magnetic held 

We assume that the magnetic energy density is much 
smaller than the matter density and pressure so that the 
magnetic-field effects can be treated as perturbations on 
a spherical non-magnetized star. The background metric 
is then given by 

ds 2 = -e 2l, dt 2 + e 2X dr 2 + r 2 {d9 2 + sin 2 ddip 2 ) , (36) 

where v and A are functions of r [6f| . The function 7 for 
the spherical stars is written as 



7 = e 



(37) 



The equilibrium state of a star is described by the set of 
the following TOV equations (6l| : 



dm , o ,„ . 
— — = 4irr*p(l + e) 
dr 



dP 

dr 



e 2X ph- 



4nPr 3 



(38) 
(39) 



dv _ 1 dP 
dr ph dr ' 

where m is defined by 

m= -(l- e - 2A ). 



(40) 



(41) 



For the unperturbed spherical stars, the gravitational 
mass M, the total baryon rest mass M*, and the internal 
thermal energy Ei nt are given by 



M = m(R) 



R 



M* = 4tt / pe A r 2 dr , 
^0 

f R 

E int = 4?r / pee x r 2 dr , 
Jo 



(42) 
(43) 

(44) 



where R denotes the circumferential radius of the star 
defined by P(R) = 0. The gravitational potential energy 
W for the unperturbed stars is defined by 

\W\ =M*+E int -M. (45) 



B. Magnetic held around a spherical star 

A profile of magnetic fields is determined by specifying 
the functional forms of f^F, T, and p, which are the arbi- 
trary functions of Following Ioka and Sasaki (23|, we 
assume that these three functions as well as the integral 
constant, C, are given by 



f = L*, 

C = Co + C2 , 



(46) 
(47) 
(48) 
(49) 



where $72, Co, C%, C2, and L are constants. Because we 
consider weak magnetic fields around a spherical star, it 
is useful to introduce a smallness parameter 77 for which 
^ = 0(rj). For the constants appearing in Equations 
(l46l) - (|49|) . we further assume that 

2 = 0(1) , Co = 0(1) , C 1 = 0(77) , 

C 2 = 0(r i 2 ),L = 0(l), (50) 

Co is determined in the background equation for Equa- 
tion (|24p. Then, up to the first order in 77, F^ v is written 
as 



F„ 



( 



-e^icscfl* 

y -a r * 



\-v 



\ 

Lcsctf* a r * 

dg^ 

-d e ^ J 



and Equation (f25j) becomes 



J 3 = phC x + 



Anr 2 sin 9 



(51) 
(52) 
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Note that the contribution of ^2 is neglected because it 
is a higher-order quantity. In terms of J 3 is given by 



J 3 = ^-W a F 3a 
1 



4 we v+\ r 2 sm 2 Q 



u+X 







sm f 



d 
dr 
1 d 



^ [ e v - x — * 
or 



89 \ sin 9 89 



(53) 



From Equations (|52|) and (|53j) . we obtain the master 
equation for the flux function ^> (the GS equation), 



,A-i. 







s 



1 a 



+47rr 2 sin 2 6phe 2X Ci + e 2(x -^L 2 ^ = . 



(54) 



Because it is the azimuthal component of the vector po- 
tential, the flux function >F may be expanded by the vec- 
tor harmonics with the axial parity as 



AndY.^M^sme—Ptie), (55) 
1=1 



where Pi is the Legendre polynomial (see, e.g., Ref. 62]). 
Substituting Equation (1551) into the GS equation fl5D) 
yields 



dr 2 



d(v-X) 2(l + l)\dipi 
dr r I dr 



XX-u) £ 2 + 



1(1 + 1) 



(1 



, d(v-\)l + l 
dr r 



-phe 2X 5l 



0. 



(56) 



Regular solutions of Equation (|56|) near the center of 
the star can be written as 

ipi = a + a 2 r 2 H , (57) 

where ao and a 2 are constants with 02 given by 

e-**L 2 -2(l + l)[(l + i)\2-V2] 



a.2 



2(21 + 3) 



■ «o 



f 

To 



p h 5} 



(58) 



Here, constants vq, v 2 , A2, po and ho are defined in the 
power series expansion of the background quantities near 
r = as follows: 



v = u + v 2 r + ■ 
A = X 2 r 2 + ■■■ , 
h = h + h 2 r 2 + 
P = Po + P2r 2 + 



(59) 
(60) 
(61) 
(62) 



Following Ref. [23j, we focus on magnetized stars 
whose exterior is vacuum and whose surface has no mag- 
netic field. At the surface of the star, then, we need to 



require two boundary conditions for the flux function, 
given by 



/ n n 

ipi = o, ~r = ■ 

dr 



at r = R . 



(63) 



where R is the radius of the unperturbed star. For the 
I 7^ 1 cases, Equation (|5"6")l becomes a homogeneous equa- 
tion. In general, then, the three boundary conditions, the 
regularity condition at the center of the star and the sur- 
face boundary conditions given by Equation (|63[) . cannot 
be satisfied simultaneously. In other words, we have to 
require ipi = for I 7^ 1. For the I — 1 case, due to the 
two boundary conditions at the stellar surface, given in 
Equation (|Q3[) . the GS equation becomes an eigenvalue 
equation with respect to the two parameters ao and L. 
The other parameter C\ can be assigned freely and deter- 
mines the strength of the magnetic fields. The remaining 
constant C 2 is related to the pressure perturbation as 
discussed below. 

The r and 9 components of the vector potential 
may be obtained straightforwardly. If we set Ag = by 
using a gauge degree of freedom, Fi 2 is given by 



-d e A r 

,A-i/ r _2 
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= A7rC 1 e x - u Lr^ 1 (r)—P 1 (9). (64) 

Requiring the regularity on the symmetry axis, we have 

A r = -4ttCi e x ~ u Lr Vi (r)P x (9) . (65) 

If another gauge condition is required, we may make a 
gauge transformation 



A h 



— d^,{f(r)Pi(0)} 



(66) 



to obtain the vector potential that satisfies a required 
gauge condition. 



C. Stellar deformation due to the magnetic fields 

As discussed before, the flux function in the present 
situation is given by 

d 



* = 47rCiT-tyi(r)sin0— P 1 (6) 
= 4 7 rC 1 * 1 (r)(|p 2 (0)-! 



(67) 



Thus, the energy-momentum tensor associated with the 
electromagnetic fields T( em )^, defined by 



, (68) 



induces a deviation of the order of 0(rj 2 ) from the back- 
ground spherical matter distribution. The line element 
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is then perturbed as follows: 



W 2 = -Le x - V {^>i) 2 , 



(76) 



ds 2 = -e 2lJ I 1 + 2 (4 7 rC 1 ) 2 F i (r)P i (e) \ dt 2 



i=0,2 
2A 



+ e 2A I 1 + 2— £ (4 7 rCi) 2 M l (r)P J (0) I dr 2 

[ r i=0,2 J 

+ r 2 {1 + 2{A^C l ) 2 K 2 {r)P 2 {6)} (d9 2 + sin 2 #d<^ 2 ) 
+ 2(4Trd) 2 W 2 (r) sin 6d e P 2 (Q)drdip, (69) 

where M i; K u I,, V i: Wi = 0(1) because Ci = 0(?/). 
Here, we employ the Regge- Wheeler gauge. In this per- 
turbed spacetime, the function 7 in Equation (|18j) is given 
by 

7 = Hl-£ (4 7 rC 1 ) 2 J ff l (r)F l (0) 1 (70) 

because Ju" 4 = and (5m 3 is the second-order quantity 
(see Equations (IT71) and (|4l))) ). Here, we have omitted 
the terms higher than 0(i] 2 ). From Equation (f24)) . the 
pressure perturbation (5P is, in terms of the metric and 
the flux functions, written as 



SP = (4^d) 2 SPi(r)Pi{ff) 
1=0,2 

= ph \ Ci* - C 2 - (47rC 1 ) 2 P l (r)P l (0) 



i=0,2 



Thus, we have 

SP = ph 
SP 2 = ph 



- — *i(r) - H (r) - C 2 



— *i(r) - P 2 (r) 

67T 



(71) 

(72) 
(73) 



where C 2 ee C 2 / (And) 2 . 

We then obtain a set of the metric perturbation equa- 
tions as follows: 



^ = 4.r> ( d{P t Pe) \ ( S 4) + V 2A 
dr f \ dP I \ ph I 3 \ dr 



1 

"3^ 



(2 + L 2 r 2 e- 2v ) (*i) : 



(74) 



dHo 
dr 



4irre 2X ph 



SPo 
ph 



g2A 

g2A 
'3^3 



1 /d$l 

3r \ 



2 + L 2 r 2 e - 2 ^ (*i) 2 , 



(75) 



er M2 2e~ 2A /d*i\ 2e 



3 V dr 



-2// 



L 2 (*0 2 , (77) 



U^ + d^ )+ d,dK 1 _ 2 e^ 
r \ dr dr / dr dr r z 

e 2A e 2A / diA 

2^ 2 - -3- 1 + 2r— M 2 

r z r J \ ar 



4ne 2X 5P 2 



1 



3r 2 V dr 



„2A 



-l I (4 + LVe^)(* 1 ) 2 , 



dP 2 di^ 1 / di/ 
dr dr r I dr 



(78) 



e 2A / diA 4 T d*i 

{ 1 + r d^J Al2 = — ^i-JT- 



3r 2 dr 



For the perturbation with / = 0, it is convenient to 
obtain Mq and 5Po/(ph) first. The equation for deter- 
mining 5Pq I (ph) is derived from Equations (|72")l and (|75|) 

as 



d_ /SPo 
dr \ ph 



1 d*! dP 



677 dr dr 

i«_wv.(^ 

67r dr \ ph 

_!_ f 1 2r— ^ M - — ( — 

r 2 \ dr / ° 3r \ dr 



.,2 A 



3r 3 



(-2 + L 2 r 2 e' 2l, )(^ 1 ) 2 . (80) 



Following Ref. 23], a new dependent variable Y 2 is, for 
the perturbation with 1 — 2, introduced by 



Y 2 ee H 2 + K 2 



-2A 



(81) 



6 



— +-^1— + — (*i) 2 
dr I r dr r A 



which facilitates the numerical computation. Then, two 
independent variables Y 2 and H 2 may be determined by 



dY 2 dv rph ( d*i 

IF = - 2 ^-^ 2 * 1 + r lF 



2 _ 2 „ r2 d^ /T 

3 6 dr" ( l} 



2,-,2-„ T , ,2 , g _ 2A ^^j (82) 



1 

3^" 



-2 + 2e" 2A <! 1 + r 



+e- 2 »L 2 r 2 } *! 



dr \ dr 
d(i/ + A) 

dr 
d*i 
dr ' 
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dr dr 



-,2A 



Aire 2X hp S H 2 



d*i 

dr 



+ 




(83) 
3r 2 dr 1 dr 



2? 



dr / 



Once ^i, 5Po/{ph), Mo, Y 2 , and i?2 are obtained, the 
other perturbation quantities may be calculated alge- 
braically through Equations ([72]). (1731) . ([75]). ([77]) . and 

Near the center of the star, the physically acceptable 
solutions may be expanded in the power series of r as 

SP /(ph) = h 00 + h 02 r 2 + • • • , (84) 

M = r 3 (™ o + mo 2 r 2 ---), (85) 

Y 2 = r 4 (y 20 + y 22 r 2 + •■•), (86) 

#2 = r 2 (h 20 + h 22 r 2 + •••)■ (87) 

Here, we have two options for determining a value of hoo- 
One is to set hoo = 0, which corresponds to considering 
sequences of the magnetized stars characterized by the 
fixed central density. The other is to use hoo to keep the 
total baryon mass M* constant for the magnetized stars. 
This corresponds to considering sequences of the magne- 
tized stars characterized by the fixed total baryon mass. 
Following Ref. [23j], we choose the latter option, i.e., on 
the constant baryon mass sequences of the magnetized 
stars. 

Outside the star, the master equations become 



dM 
dr 



0. 



dHo 
dr 



(r - 2M) 2 ' 



W 2 = . 



H 2 + 



M 2 
r-2M 



= 0. 



2M 



dY 2 = 

dr _ r(r - 2M) 



H 2 , 



dH 2 
dr 



2 \r + r{r - 2M)\ H2 M* 2 



(89) 
(90) 
(91) 

(92) 
(93) 



For the I = perturbation, the vacuum solutions are 
given by 



Mo = M (R) = const. , H Q 



Mp{R) 
r -2M 



(94) 



As for the I — 2 perturbation, manipulating Equations 
and flU]) yields 



-(y + l)(y-l) 
G 



dy 2 
4 



2y- 



dHi 



(y + 1)0,-1) 



dy 
H 2 



(95) 



where y = r/M — 1. This is the associated Legendre 
equation. Since lim H 2 = for physically acceptable 

solutions, we obtain 



H 2 = DQ 2 {y), 



(96) 



where Q™ and D are the associated Legendre function 
of the second kind and a constant, respectively. With a 
recurrence relation for Q™, 



dQ\ n 
dy 

we have 



my 



Yi = - 



2D 



--Ql(y)- 



(97) 



(98) 



At the surface of the star, the outer solutions, given by 
(IM1) . and (|98|) . are matched to the inner solutions 
integrated from the center of the star with the boundary 
conditions ~ ([57|>. 



D. Global quantities characterizing magnetized 
stars 

As mentioned before, the global physical quanti- 
ties (|3"3"]) - (|3"5T) are useful for exploration of the magne- 
tized star. Perturbations due to the magnetic effects in 
the gravitational mass, the total baryon rest mass M*, 
and the internal thermal energy E- m t are, respectively, 
given by 



AM = (47rC*i) 2 M (i?), 
AM* = 47r(47rCi) 2 



pe x r 2 



AE int = AnUM 



dP 



2,\ 



M 



pee x r 2 



dln(pe) 
dP 



SP 



s 2A M 



(99) 



dr , (100) 



dr (101) 



As already mentioned, we study the sequences of equilib- 
rium states of the magnetized star characterized by the 
fixed total baryon mass. Thus, the condition of AM* = 
is employed for determining values of hoo ha Equation 
([84]) . which is related to the perturbations in the central 
density of the star, Ap c , through the relation 



A Pc = (AnC^ph^ 



r=0 



l 00 • 



(102) 
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The electromagnetic energy £em is decomposed as 

^ (t) (103) 



EM 



EM ' 



where and fi^ are the poloidal and toroidal 

magnetic-field energies, respectively, given by 



E 



(p) _ 


16tT 2 C 2 


EM — 




f 




JO 




,(t) _ 


Wtt 2 C 2 


EM — 


3 



,A„2 .1.2 



'r 4 i>l dr . 



dr ,(104) 



(105) 



Multipolc moments are also global and physical quanti- 
ties characterizing the equilibrium star. The constant D 
appearing in the outer solutions, Equations (|96| and (1M1) . 
is related to the mass quadrupole moment AQ, defined 
by AQ = 8M 3 D/5 (see, e.g., Refs. HHH). 

Characteristic quantities for the stellar deformation 
due to magnetic stress also feature the magnetized 
stars. The surface of the star is defined by P(r) + 
(47rC*i) 2 (<5P (r) + 5P 2 (r)P 2 (9)) = 0. Thus, the radial 
displacement of the fluid elements on the stellar surface, 
Ar, is given by 



Ar = (Ar) + (Ar) 2 P 2 (9) , 



dr 



= -(AndfiSPoiR) + SP 2 (R)P 2 (e)) — (R\106) 

Here, (Ar)o may be interpreted as an average change in 
the radius of the star induced by the magnetic effects. 
The degree of the quadrupole surface deformation due 
to the magnetic stress is well described by the ellipticity, 
given by 



--(4ttC 1 ) 2 



(Ar) 2 
R 



K 2 (R) 



(107) 



where e* is defined as a relative difference between the 
equatorial and polar circumference radii of the star [23j j . 
Thus, e* < (e* > 0) means that the star is prolate 
(oblate). 

Another physically important quantity of magnetized 
objects is the total magnetic helicity H, which is con- 
served in ideal magnetohydrodynamics, and is defined 
by 



u 



(108) 



where H° is the time component of the magnetic helicity 
four-current H^, defined by 



H» = ~e!°"*f > A v F a f ) . 



(109) 



Taking the covariant derivative of Equation (|109j) yields 

1 



2 



(110) 



0, i.e., F*» v = B^u v 



Thus, V^iP 1 = if F^u" 
B^u^, and we confirm that the magnetic helicity H is a 
conserved quantity in ideal magnetohydrodynamics. For 
the present models, the total magnetic helicity is explic- 
itly written as 



u 



16tt 



(47tCi) 2 L 



(111) 



where the surface boundary condition (1631) has been 
used. The dimensionless magnetic helicity, defined by 
TIm = T-L/M 2 , is used when its numerical value is shown. 
The magnetic helicity is a measure of the net twist of a 
magnetic-field configuration. Thus, the magnetic helicity 
vanishes for purely poloidal fields and for purely toroidal 
fields. Some experiments and numerical computations 
show an interesting fact that the total magnetic helicity 
is likely to be conserved even when the resistivity cannot 
be ignored [10, [63[ . If this fact is retained for the neu- 
tron star formation process, the total magnetic helicity 
has to be approximately conserved during its formation 
process. 



V. NUMERICAL RESULTS 

In this section, we present some numerical examples 
of stably stratified stars composed of mixed poloidal- 
toroidal magnetic fields. As one-parameter equations of 
state, we employ the following one, 



P 



K p 



1+1 
1 P 



r-ip' 



(112) 
(113) 



where k and n are the polytropic constant and index, 
respectively, and Y denotes the adiabatic index, which is 
defined by Y — (din P/ din p)s- k, n, and Y are constants 
and may be specified independently for the construction 
of equilibrium stars, i.e. Y is not 1 + 1/n in general. We 
define a general relativistic Schwarzschild discriminant A 
for non-magnetized spherical stars by 



.4 



1 dp* 11 dP 
ph dr Y P dr 



(114) 



where p* is the total energy density, defined by p* = 
p + pe. Note that there is no unique definition for the 
general relativistic Schwarzschild discriminant. However, 
in the present context, only its sign matters. For different 
definitions of it, see, e.g., Refs. |64Tj66| . Following Ipser 
and Lindblom [66j, we employ a definition of the Brunt- 
Vaisala frequency N = \J—gA with 
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2(«/-A) 



l_dP 

ph dr 



(115) 



If there is a region of A > 0, the star has an unsta- 
bly stratified region and becomes convectively unstable 
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there. For the standard stars whose density profile ev- 
erywhere satisfies dp/dr < 0, the condition of the stable 
stratification for the equations of state (|112j) and ()113[) 
is given by 



Note that for the isentropic case, defined by T = (n + 
l)/n, the star is not stratified (marginally stable against 
convection), as already mentioned. 

Hereafter, we consider the n — 1 case for simplicity, 
while we employ T = 2 and 2.1: The T = 2 models 
are non-stratified ones, whose results are compared with 
those given in Ref. [23] to check our numerical results. 
The choice of T — 2.1 comes form the following reason: 
As argued by Reisenegger and Goldreich [5l[ , the Brunt - 
Vaisala frequency N from the proton-neutron composi- 
tion gradient inside neutron stars is approximated by 

where x denotes the ratio of the number densities of pro- 
tons to neutrons. From Equation (|117l) . it is found that 
the r = 2.1 models have Brunt-Vaisala frequencies at the 
stellar center similar to those of Reisenegger and Goldrc- 
ich's models with x ~ 0.1, which is a reasonable value 
for sufficiently cold neutron stars. Note that the origins 
of the buoyancy in our models and normal neutron stars 
are different. In our models, the buoyancy results from 
the entropy gradient not the composition gradient. 

First, we describe properties of the unperturbed stars. 
In the Newtonian framework, the profiles of a star like 
a density distribution are independent of values of T be- 
cause Newtonian gravity is determined only by the rest 
mass density and the pressure is assumed to be a func- 
tion of the rest mass density only. Thus, we may calculate 
thermodynamical structures of the star like internal en- 
ergy for any value of T after determining the structure 
of a star. In general relativity, by contrast, the profiles 
of the star do depend on the value of T because P and 
e are a source of gravity. Thus, we have to recalculate 
stellar profiles in general relativity whenever a value of T 
is changed. 

In Figure [TJ the gravitational mass M and the baryon 
mass M* are plotted as functions of the central density 
of the star go = p( r = 0). Throughout this paper, we 
use units of k = 1 when showing numerical results. This 
figure shows that values of M and M* for the T = 2.1 
models are larger than those for the r = 2 models for 
the same central density. Figure [2] plots relative differ- 
ences between the T = 2.1 and the T = 2 models in 
the gravitational mass M, the baryon mass M* , and the 
stellar radius R, as functions of the central density qo- 
This shows that the radius of the star increases with in- 
creasing r if one keeps the central density constant. The 
reason is that an increase in T reduces the total internal 
energy of the star and leads to a decrease in an effective 



TABLE I: Global and physical quantities for the background 
stars in units of k, — 1. 



r 


M/R 


qo 


M 


M* 


E int /\W\ 


2.000 


0.1000 


0.07027 


0.1062 


0.1118 


0.4078 




0.2000 


0.25582 


0.1623 


0.1780 


0.5633 


2.100 


0.1000 


0.06983 


0.1066 


0.1124 


0.3706 




0.2000 


0.24893 


0.1643 


0.1820 


0.5104 



TABLE II: Dimensionless low-order eigenvalues Li in units of 

iT\ 

T Mode order M/R = 0.1000 M/R = 0.2000 



RL 1 


5.792 


3.907 


RL 2 


8.315 


5.601 


RL 3 


10.80 


7.257 


RL 4 


13.26 


8.903 


RL 5 


15.71 


10.54 


RL 6 


18.16 


12.18 


RLi 


5.787 


3.909 


RL 2 


8.316 


5.619 


RL 3 


10.80 


7.283 


RL4, 


13.26 


8.937 


RL 5 


15.71 


10.59 


RLo 


18.16 


12.23 



gravitational attraction force. Figure [2] also shows that 
SM* > SM > 8R for all the values of qo calculated in 
the present study and that a typical value of the relative 
difference in the baryon mass for a standard neutron star 
model with go ~ 0.2 is within several percents. Thus, 
we may conclude that the r — 2.1 model has proper- 
ties quite similar to those of the T = 2 model, except 
for the stability against convection: We emphasize again 
that the stratification condition of the T = 2.1 model is 
absolutely different from that of the T = 2 model. Fig- 
ures [T] and [5] show that the gravitational mass, baryon 
mass, and radius of the stars tend to be independent of 
r as qo — > 0, i.e., in the Newtonian limit. Henceforth, 
we will focus on the compact models with M/R = 0.1 
and 0.2. The models having M/R = 0.2 will be a rea- 
sonable model of a neutron star. In Table U some global 
and physical quantities for the r = 2 and 2.1 models are 
tabulated. 

Next we explore the magnetic-field profile around the 
background spherical stars and its effects on the stellar 
structures. As already mentioned in the previous sec- 
tion, the magnetic-field profile is determined by solving 
an eigenvalue equation with respect to the eigenvalue L, 
which is related to the toroidal magnetic-field strength 
through L = F\2 sin 6>e ly ~ A, I , ~ 1 . Then, a discrete se- 
quence of L is allowed for the magnetized stars satisfying 
the boundary conditions ([57)) and ([55]) . Table [H] lists a 
sequence of the dimensionless eigenvalues L in units of 
R^ 1 for the first six eigensolutions, where Li means the 
i-th eigenvalue L satisfying < L\ < L 2 < L 3 < ■ ■ ■. For 
the models with T — 2, the values of L in units of 
are compared with those given in Table I of Ref. [23[ and 
we confirm that our results are in excellent agreement 
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with theirs. 

Here, we should remark the following point. In the 
study of Ref. 23], the meridional flow is in general 
present. However, the meridional flow may be absent in 
the solutions of Ref. (23[. This no meridional- flow limit 
is given by the limit |C| — > oo with L kept to be con- 
stant, in the notation of Ref. [23T ]. Even if we take this 
limit, the magnetic field remains unchanged (see Equa- 
tions (83) - (86) of Ref. [23|) although the meridional flow 
vanishes (see Equations (88) and (89) of Ref. (H). By 
using a relation he v — const., which is satisfied for non- 
magnetized "isentropic" spherical stars, we also show for 
isentropic stars with no meridional flow that basic equa- 
tions of Ref. |H, (75), (94), (95), and (100) - (103), 
become equivalent to our basic equations, (|5rJ|) . and ([74]) 

(|79p (but note some differences in notation and def- 
inition of physical quantities, e.g., their L corresponds 
to our L). Here, it should be emphasized that the pa- 
rameter C associated with the strength of the meridional 
flow does not appear in the equations of Ref. 23], (75), 
(94), (95), and (100) - (103), which fully determine the 
structure of magnetized stars with no meridional flow. 

TablcHTlshows that the dimcnsionless eigenvalues RxL 
for the r = 2.1 model are approximately equal to those 
for the r = 2 model. The same result is found in the flux 
function if one regards ^(r, 9) as a function oir/R not 
as r. These facts imply that slight changes in T do not 
affect values of R x L and ^(r/R, 9) distributions. 

Figures [3] and 0] display the profiles of magnetic fields; 
eigensolutions of "J and Fi 2 with L = L\ - L 6 for the 
spherical star with M /R = 0.2 and r = 2.1 (see the cor- 
responding eigenvalues in Table HI)) . Figure [3] shows how 
lines of the magnetic force on the meridional cross sec- 
tion behave, because an equi-^ line corresponds to a line 
of the magnetic force. These figures suggest that higher- 
order eigensolutions have more non-uniform structures of 
the magnetic fields. Figure [3] also shows that there is a 
negative region of for the models with L 2 , L 4 , and L 6 , 
while '5 is everywhere non-negative for the models with 
L = Li, L3, and L5. 

For analyzing properties of magnetic-field profiles, it is 
helpful to introduce an orthonormal tetrad component of 
the magnetic field, -B M , given by 

£?( t ) = , £?( r ) = — SttCi^x cos 9 , 
B {0) =47rC*ie- A (r^i + 2?/>i)sm(9, (118) 
-8(0) = 47rCie _iy Lripi sin 9 . 

Then, we can define that B c = \B(r = 0)| = 87^1^1(0). 
As mentioned before, profiles of for the T = 2.1 
models are very similar to those for the T = 2 models. 
Their difference cannot be visible if are plotted as 
functions of r/R in the same figure even though no figure 
is given in this paper. To show the dimensionless total 
magnetic helicity of the magnetized stars, Umi we follow 
Ioka and Sasaki [23j and use a dimensionless parameter 



that represents magnetic-field strength, defined by 



which is as large as the ratio of the magnetic energy to 
the gravitational energy. 

Table Hill lists global physical quantities characterizing 
the magnetized stars with mixed poloidal-toroidal fields; 
the changes in the central density Ap c , the gravitational 
mass AM, the quadrupole moment AQ, the mean ra- 
dius (Ar)o, the ellipticity e* , the magnetic helicity 'H.m-, 
the magnetic energy -Eem, and the ratio of the poloidal 
magnetic energy .E^m ^° ^ ne total magnetic energy Eem- 
In this table, the results for the first six eigensolutions 
are shown and all the quantities are normalized to be 
non-dimensional, as given in the first row. By comparing 
the results shown in Table IIIII with those in Table 2 of 
Ref. [23j], we can again check the reliability of our results; 
ours are in agreement with theirs for the T = 2 models 
in an acceptable level. 

Although there are slight numerical differences be- 
tween them, there is no qualitative difference between 
the results of the r — 2.1 and T = 2 models (see Ta- 
ble [TTTJ) . Thus, basic properties of the magnetized star 
described in Ref. [23[ hold for the present r = 2.1 mod- 
els, although the convective stability is different and also 
any meridional flow is absent in the present models. For 
the magnetized-star sequences, characterized by the fixed 
baryon mass and magnetic helicity, shown in Table IIIII 
the total gravitational mass and the total magnetic en- 
ergy increase with the mode number i (or the eigenvalue 
L). This property is reasonable due to the following rea- 
son: As found in Figures [3] and 01 higher-order eigensolu- 
tions, characterized by a lager eigenvalue, have more non- 
uniform structures of the magnetic fields, which can usu- 
ally store larger magnetic-field energy. From Table IIIII 
we also find that all the models obtained in this study 
are toroidal-magnetic-field dominant and that the value 
of E^^/Eem decreases as the mode number i increases. 
The later property is well explained by the fact that L 
is interpreted as the ratio of the toroidal magnetic-field 
strength to the poloidal magnetic-field one, and it in- 
creases as the mode number i increases. The magnetic 
hoop stress around the symmetry axis due to the toroidal 
magnetic field tends to make the star prolate like a rub- 
ber belt fastening a waist of the star. This is consistent 
with the facts that all the stars obtained in this study 
have negative ellipticity e*, i.e., the star is prolate, and 
that the degree of the quadrupole deformation measured 
by |e*| becomes more pronounced for higher-order solu- 
tions, as shown in Table HTT1 

VI. DISCUSSION: THE STABILITY OF THE 
STARS 

Checking the stability is an important issue to be ex- 
plored after obtaining an equilibrium state of a magne- 



12 



TABLE III: Global and physical quantities characterizing the magnetized stars with mixed poloidal-toroidal fields; the changes 
in the central density Ap c , the gravitational mass AM, the quadrupole moment AQ, the mean radius (Ar)o, the ellipticity e*, 

the magnetic helicity Hm, the magnetic energy -Eem, and the ratio of the poloidal magnetic energy J?^m *° ^ ne total magnetic 
energy Eem- Here, all the quantities are normalized to be non-dimensional, as given in the first row. 



[V, M/R) 


L 


Apc/qoHM 


AM/MHm 


AQ/MR 2 Um 


(Ar) /RH M 




e* /Hm 






Hm 


/Km 




fiEM I\W\Hm 




(2,0.2) 


Li 


0.2545 


2.446 


X 


10" 


2 


-6.163 


X 


10" 


d 


-2.601 x 10 


— 2 


— 1 


505 


X 


10" 


2 


2 


556 


X 


10" 


i 


0.1903 


0.3684 




Li 


0.4063 


3.067 


X 


10" 


2 


-1.875 


X 


10" 


2 


9.007 x 10" 


2 


—4 


580 


X 


10" 


2 


4 


346 


X 


10" 


i 


0.2222 


0.2785 




L 3 


0.5307 


3.634 


X 


10" 


-2 


-3.293 


X 


10" 


2 


4.156 x 10" 


2 


—8 


043 


X 


10" 


2 


1 


153 


X 


10" 


i 


0.2606 


0.2097 




Li 


0.6736 


4.202 
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tized star, because unstable solutions lose their physical 
meaning in the sense that they are not realized in na- 
ture. For the present models, as observed in Sec. V, 
the gravitational mass and the total electromagnetic en- 
ergy increase with the mode number i of the eigensolu- 
tion, if we consider the equilibrium sequences for a given 
baryon mass and a magnetic helicity. This result is quite 
important and interesting due to the following reason. 
If the total baryon mass and magnetic helicity are con- 
served during the formation process of neutron stars, as 
discussed before, it is likely that the final state of the 
magnetic fields characterized by the arbitrary functions 
of the flux function P5|) - (|4"5)l and the surface boundary 
condition (|63|) is the lowest-order eigensolution charac- 
terized by the smallest eigenvalue L = L\ because it has 
the lowest gravitational mass and total electromagnetic 
energy among the equilibrium solutions characterized by 
the fixed baryon mass and magnetic helicity. We there- 
fore conjecture that all the high-order solutions are un- 
stable because there is an equilibrium state with energy 
lower than theirs and that the solutions associated with 
the lowest eigenvalue L = L\ can be stable if there is a 
stable solution among the present models. Another im- 
portant fact is that for the magnetic-field profile charac- 
terized by L = L±, their magnetic energy is most equally 
divided into the poloidal and the toroidal magnetic en- 



ergies among the eigensolutions. This property is con- 
sistent with the conjecture for stable magnetic config- 
urations given by the linear analyses; stable magnetized 
stars contain both poloidal and toroidal components with 
comparable magnetic energies. 

Most of the magnetized-star models constructed in 
the framework of general relativity so far were non- 
stratified, and therefore, marginally stable against con- 
vection. Those non-stratified models are in general highly 
unstable against the magnetic buoyancy in the vicinity of 
the stellar surface, in the presence of magnetic fields. For 
strongly magnetized stars like magnetars, as shown by 
Kiuchi et al. (48j, the magnetic buoyancy instability in- 
duces a convective motion near the surface of the star and 
fully destroys initially coherent magnetic fields inside the 
star. To stabilize this magnetic-buoyancy instability, the 
stratification with the strength sufficient to overcome the 
magnetic buoyancy are necessary as a stabilizing agent. 
This stabilization effect prevails in non-rotating diffusion- 
less stars wit h p urely toroidal magnetic fields as argued 
by Acheson j35|. When A^ 2 lo\ > 0, his dispersion 
relation (see Equation (3.20) of Ref. [35|]) has four so- 
lutions, uj « ±kgN(kr + kg)~z and lu « zLuilua, where 
lu and oj^ mean the oscillation and the Alfven frequen- 
cies, respectively, and k r , kg, and m denote the ver- 
tical, horizontal, and azimuthal wave numbers, respec- 
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tively. These four solutions are composed of propagating 
waves; the former is an internal gravity wave and the lat- 
ter an Alfven wave. We therefore confirm that there is no 
magnetic-buoyancy instability as long as N 2 uj\ > 0. 
In Figure [5j we plot squares of the Brunt- Valsala, fre- 
quency N 2 and the Alfven frequency ui\ for a T = 2.1 
model characterized by M/R = 0.2 and L = L\ as func- 
tions of the dimensionless radius r / R. Here, the Alfven 
frequency is evaluated on the equatorial plane and de- 
fined by uj a = y/B^B^/HAirph + B^B^r 2 ), and the 
strength of the magnetic fields is determined by the con- 
dition -Eem/IW^I = 2.5x 10 -2 , which corresponds to very 
strong magnetic fields « 10 16 G for a typical neutron star. 
This figure shows that the Brunt- Vaisala frequency, N, 
is much larger than the Alfven frequency uja in the re- 
gion of r > 0.5R for the models with L — L\, Thus, 
we can predict that this model is stable against the mag- 
netic buoyancy vicinity of the stellar surface. It should be 
noted that the physical origin of the anti-buoyancy force 
in our models is different from that in real neutron stars, 
as mentioned before, because the former comes from the 
entropy gradient and the later mainly from the composi- 
tion gradient (see, e.g., Ref. [5lj|). Magnetic fields, how- 
ever, do not care the origin of the anti-buoyancy force, 
but do whether stars are stably stratified or not. 

As pointed out by Acheson (36|, the magnetic hoop 
stress caused by the strong toroidal magnetic fields gov- 
erns dynamics of the perturbation near the center of 
the star. This fact may be confirmed from Figure [SJ 
which shows N <C uja near the center of the star. Thus, 
the stratification is not helpful near the center of the 
star in the presence of the magnetic instability. For the 
central part of the star, as mentioned before, the pres- 
ence of the poloidal magnetic fields having comparable 
strength with the toroidal ones will suppress the pinch- 
type instabilities. By solving linear perturbation equa- 
tions around magnetized star models in the framework of 
Newtonian dynamics, Lander and Jones showed that this 
suppression of the magnetic instability indeed occurs, al- 
though the presence of the poloidal magnetic fields leads 
another instability associated with themselves [46| (see, 
also, Refs. [H, Thus, the results of their numerical 
simulation lessen the possibility that the pinch-type mag- 
netic instabilities are completely removed for the stars 
containing the mixed poloidal-toroidal magnetic fields 
with comparable strength. Although Lander and Jones's 
results obviously conflict with those by Braithwaite and 
his collaborators (40H431 ]. we have so far had no definite 
conclusion to this controversy. One important difference 
between their analyses is in the treatment of the resistiv- 
ity of the matter; Lander and Jones employed the ideal 
magnetohydrodynamics approximation, whereas Braith- 
waite and his collaborators took into account the resistiv- 
ity. This might be a key to the solution of this problem. 
Toward a definite answer to this stability problem, we 
need further studies on the stability of the magnetized 
stars. 

As such a study, we plan to perform a stability anal- 



ysis of the present magnetized star models by general 
rclativistic magnetohydrodynamics simulations like those 
done by Kiuchi et al. [13, EH . The present stably strati- 
fied models characterized by L = L\ satisfy the following 
conditions; they have both poloidal and toroidal mag- 
netic fields with comparable strength to suppress the 
hoop-stress instability inside the star, and in addition, 
the fluid constructing the magnetized star is stably strat- 
ified with strength sufficient to overcome the magnetic 
buoyancy near the stellar surface. Thus, it is obvious 
that some major magnetic instabilities will be reduced in 
the present magnetized star models. This will make the 
stability problem more tractable. 

As discussed in Kiuchi et al. [U E3| , it is reasonable 
to expect that the toroidal component of the magnetic 
fields is much larger than the poloidal component in- 
side the neutron star at least soon after its birth. The 
reason for this is that the winding of poloidal magnetic 
fields caused by a rapid and differential rotation dur- 
ing the core collapse would create large toroidal fields. 
Most of general relativistic magnetized star models ob- 
tained in numerical computations so far can, however, 
have toroidal magnetic fields much weaker than poloidal 
ones [2~il - |2"6l | . Their minimum ratio of the poloidal mag- 
netic field energy to the total magnetic energy is « 0.92, 
which is much larger than those of the present models. 
(For magnetized star models in the framework of the 
Newtonian dynamics, see, e.g., Ref. [67j). Their mag- 
netic fields are composed of the mixed poloidal-toroidal 
twisted torus magnetic fields inside the star and nearly 
dipolar magnetic fields outside the star. They may look 
quite plausible for the magnetosphere of neutron stars. 
However, weak toroidal magnetic-field strength seems to 
be unlikely for the strongly magnetized neutron stars. 
Although the magnetic field vanishes outside the star for 
the present models, which is quite unrealistic, the present 
models would give a reasonable inside structure of the 
strongly magnetized neutron stars because of their large 
toroidal magnetic-field strength. For obtaining more re- 
alistic models composed of a stably stratified fluid, basic 
equations given in Equations (|19p - (I25|) may be em- 
ployed as far as ideal magnetohydrodynamics and the 
barotropic equations of state are employed. 



VII. SUMMARY 

We constructed the magnetized stars composed of a 
stratified fluid in the framework of general relativity. By 
assuming ideal magnetohydrodynamics and employing a 
barotropic equation of state, we first derive basic equa- 
tions for describing stably stratified stationary axisym- 
metric stars containing both poloidal and toroidal mag- 
netic fields. As sample models, the magnetized star con- 
sidered by Ioka and Sasaki [23| are modified to the ones 
stably stratified. The resulting models have both poloidal 
and toroidal magnetic fields with comparable strength. 
The magnetized stars newly constructed in this study 
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are believed to be more stable than the existing relativis- 
tic models because they have both poloidal and toroidal 
magnetic fields with comparable strength, and the mag- 
netic buoyancy instability near the surface of the star, 
which can be stabilized by the stratification, are sup- 
pressed. 
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FIG. 1: Gravitational mass M and baryon mass M* , given as functions of the central density qo- 




FIG. 2: Relative differences of the gravitational mass M, the baryon mass M* , and the radius 7?, given as a function of the 
central density qo- Here, the relative difference of the physical quantity Q[T] is defined by SQ = 2(Q[2.1] — Q[2])/(Q[2.1] + Q[2]). 
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FIG. 3: Equi-\& contours on the meridional cross section for the models characterized by L — L\, L2, £3, L4, L5, and Lq. 
Here, z and vj are defined by z = rcosO and w = rsmO, respectively. The thick quarter circle shows the surface of the star, on 
which \f = is, by the boundary condition, required. The black and white regions, respectively, correspond to the maximum 
and minimum values of \&. Since ^ vanishes outside the star, it can be seen that there is a negative region of * for the models 
with Z/2, L4, and Le, while ^ is always positive for the models with L — Li, L3, and L5. 




FIG. 4: Equi-Fi2 contours on the meridional cross section for the models characterized by L = L\, L2, L3, L4, L5, and Lq. 
Here, z and w are defined by z = r cos 6 and zu — r sin 0, respectively. The thick quarter circle shows the surface of the star, on 
which F12 = is, by the boundary condition, required. The black and white regions, respectively, correspond to the maximum 
and minimum values of F12 . Since F12 vanishes outside the star, it can be seen that there is a negative region of _Fi2 for the 
models with L2, L4, and Lq, while F12 is always positive for the models with L = Li, L3, and L5. 
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FIG. 5: Squares of the Brunt -Vaisala frequency N 2 and the Alfven frequency given as functions of the dimensionless 
radius r/R. Here, the Alfven frequency is defined by lua = \J B^B^/^irph + B^B^)r 2 ) and evaluated on the equatorial plane 
(8 = 7r/2), and the strength of the magnetic fields are determined by the condition -Eem/IW = 2.5 x 1CP 2 . 



